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\ Abstract 

I— I \ We report basic far-infrared (FIR) properties of eight blue compact dwarf galaxies (BCDs) observed 

i-C ■ by AKARI. We measure the fluxes at the four FIS bands (wavelengths of 65 yum, 90 /.im, 140 /.im, and 

160 /im). Based on these fluxes, we estimate basic quantities about dust: dust temperature, dust mass, 
Q . and total FIR luminosity. We find that the typical dust temperature of the BCD sample is systematically 

| H ' higher than that of normal spiral galaxies, although there is a large variety. The interstellar radiation 

^ , field estimated from the dust temperature ranges up to 100 times of the Galactic value. This confirms 

the concentrated star-forming activity in BCDs. The star formation rate can be evaluated from the FIR 
luminosity as 0.01-0.5 Mq yr~^. Combining this quantity with gas mass taken from the literature, we 
estimate the gas consumption timescales (gas mass divided by the star formation rate), which prove to 
\ span a wide range from 1 Gyr to 100 Gyr. A natural interpretation of this large variety can be provided by 

intermittent star formation activity. We finally show the relation between dust-to-gas ratio and metallicity 
(we utilize our estimate of dust mass, and take other necessary quantities from the literature). There is a 
positive correlation between dust-to-gas ratio and metallicity as expected from chemical evolution models. 
Key words: galaxies: dwarf — galaxies: ISM — infrared: galaxies — ISM: dust, extinction 
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2002). 

Dust grains absorb stellar ultraviolet (UV)-optical light The origin of the cosmic dust is important in itself, 
and reprocess it into far-infrared (FIR), thereby affect- and the formation and destruction of grains are tightly 
ing the energetics of interstellar medium. In particular, related to the star formation activity of galaxies. The con- 



■ Hirashita & Ferrara (2002) show that the dust produced densation of heavy elements is an important process for 
by supernovae efficiently absorbs UV heating photons and the formation of grains. In particular, in young galaxies, 
helps gas cooling and H2 production even in primeval dust is predominantly condensed in Type II supernovae 
galaxies (typical age is < 10^ yr). Dust grains also af- (SNe II), whose progenitors have short lifetimes (Dwek & 
feet the escape fraction of UV photons out of galaxies. Scalo 1980; Kozasa, Hasegawa, & Nomoto 1989; Morgan 
Observationally, this means that there is always a correc- & Edmunds 2002; Nozawa et al. 2003). However, how 
tion factor for dust extinction when one estimates star much dust forms in a SN II is not known. There is a 
formation rate (SFR) from UV luminosity (Buat & Xu large difference among observationally derived dust mass 
1996). Indeed, a large correction factor for dust extinc- (Moseley et al. 1989; Dunne et al. 2003; Hines et al. 
tion is required for the cosmic star formation rate traced 2004; Sugerman et al. 2006; Meikle et al. 2007; Ercolano, 
with UV (Steidel et al. 1999; Takeuchi, Buat, & Burgarella Barlow, & Sugerman 2007; Rho et al. 2008). Thus, it is 
2005; Le Floc'h et al. 2005; Perez- Gonzarez et al. 2005). important to constrain the dust mass formed in SNe II by 
Also theoretically, how dust accumulates in the history quantilying dust mass in young galaxies, 
of the Universe is important in the context of the cosmic A large amount of dust has been suggested to exist at 
reionization, because dust absorption affects the escape high redshift (z) (e.g., Bertoldi et al. 2003; Priddey et al. 

2003). The estimate of dust mass in high- 2; galaxies in- 
deed constrains how much dust forms in SNe II (Dwek, 
Galliano, & Jones 2007). The rest-frame UV extinction 
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curves also suggest that dust enrichment by SNe II is in- 
deed occurring at high-z (Maiolino et al. 2004; Bianchi 
& Schneider 2007; Hirashita et al. 2008). However, it is 
not easy to explore the first dust enrichment in primeval 
galaxies at high z (z ^ 5) with present observational facili- 
ties. Therefore, nearby templates of primeval galaxies are 
useful to test galaxy formation scenarios. The best can- 
didates for such a template are metal-poor (~ 1/10-1/3 
solar) blue compact dwarf galaxies (BCDs), since they are 
at the early stage of chemical evolution and their typical 
stellar age is young (Sargent & Searle 1970). Indeed BCDs 
harbor appreciable quantities of dust (Thuan, Sauvage, 
Madden 1999; Plante & Sauvage 2002; Hunt et al. 2006), 
and can be used to investigate the dust properties in chem- 
ically unevolved galaxies (Takeuchi et al. 2005). 

The dust content in a galaxy is often estimated by mea- 
suring the dust emission in FIR. Because of the large all- 
sky sample, the Infrared Astronomical Satellite (IRAS) 
60 /im and 100 /im bands are most frequently used. 
However, Shibai, Okumura, & Onaka (1999), from the 
study of FIR spectral energy distribution (SED) of the 
Milky Way, show that the IRAS bands fail to detect a 
significant fraction of emission from large grains which 
contribute to wavelengths (A) longer than 100 //m (see 
also Okumura et al. 2002). Moreover, with the size distri- 
bution derived for the Milky Way dust grains, the largest 
contribution to the dust content comes from large grains 
(a ~ 0.1 /^m, where a is the grain radius) (Mathis, Rumpl, 
& Nordsieck 1977). Large grains also have an importance 
that their temperature reflects the stellar radiation field 
of interstellar medium through the equilibrium between 
the heating from stellar radiation and the radiative cool- 
ing in FIR (Draine & Anderson 1985). In summary, FIR 
observation of emission from large grains at A > 100 /im 
is important to quantify the dust content and the stellar 
radiation field. 

Indeed, Galliano et al. (2005) show that some BCDs 
have a cold dust component contributing to A > 100 /im. 
Popescu et al. (2002) also point out that BCDs tend to 
have the coldest dust temperature by using a Virgo cluster 
sample observed by the Infrared Space Observatory (ISO), 
although the sample size is small. This is unexpected from 
the warm IRAS 60 /im/100 /im colors of BCDs (Hoffman 
et al. 1989). Engclbracht et al. (2008), by using the 70 
/im and 160 /im bands of Multiband Imaging Photometer 
for Spitzer (MIPS) on the Spitzer Space Telescope, show 
that the dust temperature tends to rise as the metallic- 
ity decreases from solar metallicity (Zq^ to ~ 1/10 Zq. 
However, there is a risk of contamination of stochasti- 
cally heated dust in the 70 /im band (Li & Draine 2001), 
which leads to overestimation of the large-grain tempera- 
ture. Thus, a statistical sample at A > 100 /im is indeed 
required to avoid possible contamination from stochasti- 
cally heated grains, and simultaneously to overcome the 
small number statistics of long-wavelength (A > 100 /im) 
observations. 

In this paper, we report observations of BCDs in FIR 
wavelengths including A > 100 /im. The data are obtained 
by AKARI (Murakami et al. 2007) Far-Infrared Surveyor 



(FIS) (Kawada et al. 2007). The continuous wavelength 
coverage (A ^ 50-180 /im) and 4 bands in the FIR provided 
by FIS enable us not only to quantify the long-wavelength 
(A > 100 /im) dust luminosity but also to separate the 
contribution from stochastically heated grains by using 
the short-wavelength (A ^ 70 /im) band. We also demon- 
strate that this advantage of FIS indeed makes it possible 
to discuss some quantities related to dust properties and 
dust enrichment. 

This paper is organized as follows. Observations and 
data reduction arc described in section 2. Then, the basic 
data are shown in section 3, and some quantities derived 
from those data are presented in section 4. In section 5 we 
discuss the results, and finally in section 6 we summarize 
this paper. 

2. Observations and Reduction 

The observations of eight BCDs are carried out by FIS 
onboard AKARI as one of the open time observing pro- 
grams (Proposal ID: BCDDE, PI: H. Hirashita) with four 
photometric bands of WIDE-S, WIDE-L, and N160, 
whose central wavelengths are 65, 90, 140, and 160 /im 
with effective band widths of AA 21.7, 37.9, 52.4, and 
35.1 /im, respectively (Kawada et al. 2007). The measured 
FWHMs of the point spread function (PSF) arc 37 ± 1", 
39 ±1", 58 ±3", and 61 ±4", respectively, for the above 
four bands. 

The sample BCDs are listed in Table 1. We primar- 
ily selected the sample with IRAS detections as listed 
in Hirashita, Tajiri, & Kamaya (2002) (originally com- 
piled by Sage et al. 1992 and Lisenfeld & Ferrara 1998). 
II Zw 40, Mrk 7, UM 439, UM 533, and II Zw 70 arc fi- 
nally observed. The limitation to the sample comes from 
the lifetime and visibility of the satellite. We also add 
Mrk 36 and Mrk 71 from Hopkins, Schultc-Ladbeck, & 
Drozdovsky (2002) since they matched the visibility re- 
quirement. Those galaxies are also detected by IRAS. 
II Zw 71 is also included in the sample, since it lies in 
the field of view of II Zw 70. This is the only galaxy 
which is not detected by IRAS. All the sample BCDs are 
observed with the FISOl scan sequence (slow-scan obser- 
vation for photometry), whose total area scanned in the 
four bands is 15' x 12'. Because of the faintness of the 
sample, we choose the slower scan speed (8 arcsec/sec) 
and the largest reset interval (2 sec). The observational 
parameters are also listed in Table 1. 

The raw data were reduced by using the FIS Slow Scan 
Tool (version 20070914). ^ The absolute calibration of 
the FIS slow-scan data is based on the COBE/DIRBE^ 
measurements of the diffuse sky emission from zodiacal 
light and interstellar cirrus averaged over areas of several 
slow-scan observations. Because the detector response is 
largely affected by a hit of high energy ionizing parti- 
cle, we use the local flat; that is, we correct the detec- 
tor sensitivities by assuming uniformity of the sky bright- 

^ http:/ /www. ir.isas.ac.jp/ASTRO-F/Obscrvation/. 
^ http:/ /cmbdata. gsfc.nasa.gov/product/cobc/. 
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Table 1. Observation log. 



Name 


R.A. (J2000.0) 


Dec. (J2000.0) 


Date 


Observing mode 


Sca.n speed 


Reset interval 


II Zw 40 


5'^55'"42?6 


3°23'31."8 


2006 Sep. 21 


FISOl 


8" s-^ 


2 s 


Mrk 7 


7'^28'"11?4 


72°34'23."4 


2006 Oct. 3 


FISOl 


8" s-i 


2 s 


Mrk 71 


7'^28'"42?8 


69°11"21."4 


2006 Oct. 4 


FISOl 


8" s-i 


2 s 


UM 439 


Ilh36'"36f8 


0°48'58."0 


2006 Dec. 16 


FISOl 


8" s-i 


2 s 


UM 533 


12h59'"58?l 


2°02'57."3 


2007 Jan. 4 


FISOl 


8" s-i 


2 s 


II Zw 70* 


14h50'"56?5 


35°34'17."6 


2007 Jan. 15 


FISOl 


8" s^i 


2 s 


II Zw 71* 


14h51'"14?4 


35°32'32."2 


2007 Jan. 15 


FISOl 


8" s-i 


2 s 


Mrk 36 


llh04'"58?5 


29°08'22."1 


2007 May 27 


FISOl 


8" s-i 


2 s 



* II Zw 70 and II Zw 71 lie in the same field of view. 



ness (Verdugo, Yamamura, & Pearson 2007).'^ By using 
some relatively bright sources, we confirmed that the er- 
rors caused by the smoothing procedures are well within 
10% for N60 (65 and WIDE-S (90 and 20% for 
WIDE-L (140 ^m). These error levels are consistent to or 
somewhat larger than the background fluctuation mainly 
caused by a hit of high energy ionizing particle, which 
implies that the background fluctuation is a major com- 
ponent in the errors. The above values are adopted for 
the errors in Table 1 except for the following two cases: 
(i) For the N60 (65 /xm) flux of Mrk 7, we adopt the root- 
mean-square (rms) of the background fluctuation, (ii) The 
uncertainty in the WIDE-S (90 fim) data of II Zw 70 is 
also large because of a bad pixel within the aperture radius 
(section 3.3). For N160 (160 /im), because of high fluc- 
tuation of the background, we adopt the 25% error after 
estimating the uncertainty in the background subtraction. 
When the flux is smaller than the rms of the background, 
we adopt 3 times the background uncertainty as an upper 
limit. 

Our sample BCDs are compact enough to be treated 
as point sources with the resolution of FIS (but see sec- 
tion 3.3 for Mrk 71). We follow the aperture photometry 
described in Verdugo, Yamamura, & Pearson (2007). We 
adopt aperture radii of 0.625 arcmin for N60 (60 fim) 
and WIDE-S (90 fim) and of 0.750 arcmin for WIDE- 
L (140 /im) and N160 (160 /im), and sky regions be- 
tween 2.25 and 3.25 arcmin for N60 (65 /im) and WIDE- 
S (90 /im) and between 3.0 and 4.0 arcmin for WIDE-I 
(140 /im) and N160 (160 fim). Then the aperture correc- 
tion factors (1.58, 1.74, 1.71, and 2.03 for N60, WIDE-S, 
WIDE-L, and N160, respectively) are multiplied to obtain 
the total flux.'' 

3 We referred to AKARI FIS Data User Manual Version 2 
available at http://www.ir.isas.ac.jp/ASTRO-F/Observation/. 
When we run the pipeline command ss_run_ss, we apply op- 
tions /local, /smooth, width_filter=80: The flat field is built 
from the observed sky, and boxcar smoothing with a filter width 
of 80 s in the time series data is applied to remove remaining 
background offsets among the pixels. 

* We also multiplied another factor called correction factor. It is 
an empirical factor to correct the discrepancy between the cali- 
bration of the diffuse emission (on which the present reduction 
is based) and that for point sources. This correction factor is 
empirically known to be constant. By using a different tool de- 
veloped by T. Suzuki, who adopt point source calibrations by 



Finally, color corrections are applied for the WIDE-L 
(140 /im) fluxes: the color correction factor is assumed to 
be 0.93 (the flux is divided by this factors), which is valid 
for the temperature range derived by using 140 /im and 
90 /im fluxes (22-53 K for /3 1 and 19-36 K for /? = 2; see 
section 4.1) within 3%. For the WIDE-S (90 /im) band, 
color correction would be required. The correction factor 
is 0.94 at 35 K and 1.06 at 50 K for /3 = 1, and 0.94 at 30 
K and 1.06 at 40 K. These temperature ranges cover the 
average of 140 /im-90 /xm temperature and 65 /im-90 /^m 
temperature (section 4.1). For II Zw 71, whose tempera- 
ture is signiflcantly lower than the above ranges, upward 
~ 10% correction would be required. However, we do not 
apply color correction to the WIDE-S (90 /im) flux, be- 
cause the emission in this band may not be described with 
a single-temperature emission and because the correction 
factor is not uniform within the sample (i.e., an iteration 
between the derived dust temperature and color correc- 
tion is required for each object). For N60 (65 /im) and 
N160 (160 /xm), we do not apply color correction, since 
the color correction change the flux only < 4% for our 
sample within the range of dust temperature derived for 
our sample in section 4.1. 

3. Results 

3.1. Flux 

In Tabic 2, we present the measured fluxes. All the 
BCDs are detected in N60 (65 fim), WIDE-S (90 /im), 
and WIDE-L (140 /im) bands except for Mrk 36 (section 
3.3). For the N160 (160 /im) band, a half of the objects 
are detected and only upper limits are determined for the 
other half. For comparison, in Table 2 we also compile 
the IRAS 60 /im and 100 /im fluxes if available. We adopt 
those data from Moshir et al. (1990) except for II Zw 40, 
Mrk 71 and II Zw 71 (section 3.3). For II Zw 40, we also 
present the data taken by Spitzer MIPS at A = 70 /im and 
160 /im (Engclbracht ct al. 2008). 

asteroids, we confirmed that the correction factors arc reason- 
able within the errors put in this paper. See also Kawada et al. 
(2007) for the absolute calibration. 
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Table 2. Measured fluxes. 



Name 



N60 
(A = 65 jum) 



WIDE-S 
(A = 90 Mm) 



WIDE-L 
(A = 140 //m) 



N160 
(A ^ 160 a™) 



IRAS IRAS 

(A = 60 ^im) (A ^ 100 ^^m) 



II Zw40* 
Mrk 7 
Mrk 71^ 
UM439 
UM533 
II Zw 70 
II Zw71 
Mrk 36 



6.9±0.7Jy 
0.83 ± 0.23 Jy 

3.7±0.4Jy 
0.42 ± 0.04 Jy 
0.54 ± 0.05 Jy 
0.89 ± 0.09 Jy 
0.32 ± 0.03 Jy 
0.33 ± 0.03 Jy 



6.6±0.7Jy 
0.84 ± 0.08 Jy 

3.1±0.3Jy 
0.40 ± 0.04 Jy 
0.61 ± 0.06 Jy 
0.80 ± 0.24 Jy 
0.54±0.05Jy 
0.26 ± 0.03 Jy 



3.7±0.7Jy 
l.l±0.2Jy 
2.2±0.4Jy 
0.31 ± 0.06 Jy 
0.54 ± 0.11 Jy 
0.75 ± 0.15 Jy 
1.3±0.26Jy 
< 0.30 Jy 



3.4±0.9Jy 
1.0±0.3Jy 
< 3.0 Jy 
<0.45 Jy 
<0.68Jy 
0.62 ±0.16 Jy 
1.3±0.3Jy 
<0.38Jy 



6.61 ±0. 7 Jy 
0.48 ± 0.04 Jy 
4.83 ± 0.04 Jy 
0.36 ± 0.05 Jy 
0.50 ± 0.05 Jy 
0.71 ± 0.05 Jy 



5.8±0.9Jy 
0.97±0.14Jy 
5.71 ± 0.03 Jy 
0.81±0.11 Jy 
0.54±0.13Jy 
1.24±0.12Jy 



0.23 ± 0.06 Jy 0.68 ± 0.08 Jy 



Note. 3 (T upper limits are shown by "<" (For the WIDE-L data of Mrk 36, see section 3.3). " — " indicates no data. IRAS 60 /^m 
and 100 /^m data are also shown for comparison in the last two columns. 

* II Zw 40 has also Spitzer MIPS data (5.58 ± 0.28 Jy at A = 70 lira and 3.14 ± 0.43 Jy at A = 160 iim; Engelbracht et al. 2008). 
t The discrepancy between our results and the IRAS data may be explained by different resolution (section 3.3). 



3.2. FIR Spectral Energy Distributions 

The data in Table 2 are displayed in Figure 1. For 
almost all the sources with available IRAS data, the 
AKARI and IRAS fluxes are consistent but there is a 
significant discrepancy between the IRAS 100 fxm flux 
(1.24 ±0.12 Jy) and the AKARI 90 (0.80 ±0.24 Jy) 
flux for II Zw 70. Note that the error bar for the 90 /im 
flux is large because of the bad pixel (section 2). However, 
we should note that the discrepancy seems significant 
even after considering the uncertainty in the interpola- 
tion. Another possible reason is that the IRAS flux is 
sometimes overestimated (Jeong et al. 2007). Source con- 
fusion could explain the discrepancy, but there is no evi- 
dence for a large sky fluctuation around II Zw 70 in our 
image. 

The N160 (160 yum) flux, when detected, is always near 
to the WIDE-L (140 ^m). Thus, the 160 /im flux docs 
not provide strongly independent information of the 140 
/im flux. Because of this and of a large fraction of non- 
detection, we do not use 160 ^m data in the following 
discussions. Nevertheless, it is worth stressing that the 
consistency between the 140 /J,m and 160 /xm fluxes con- 
firms the rehability of the N160 band. The N160 flux of 
II Zw 40 is also consistent with the Spitzer 160 /im flux, 
which also strengthen the reliability of the N160 fluxes. 

3.3. Remarks on individual objects 

II Zw 40 — The IRAS data are taken from Vader et al. 
(1993) (see also Beck et al. 2002). The Spitzer MIPS data 
are also available in Engelbracht et al. (2008). 

Galliano et al. (2005) investigate this galaxy by using 
detailed SED models. They suggest a very cold dust con- 
tributing to the emission at A ^ 300 fim. (see also sections 
4.5 and 5.2). However, we cannot constrain the dust emis- 
sion at such a long wavelength range. Hunt, Bianchi, & 
Maiolino (2005) also study the SED over a wide wave- 
length range including FIR, and state that the FIR SED 
has a peak at A ^ 60 //m, which is consistent with our 
data. This galaxy has the strongest interstellar radiation 
field, which supports high density of star formation (Beck 
et al. 2002; Vanzi et al. 2008). 



Mrk 71 — This is called "cometary" BCD (Noeske et al. 
2000) and given an NGC number NGC 2366. This is 
the nearest BCD of our sample with distance of 3.4 Mpc 
(Tolstoy et al. 1995). Indeed, an extended structure in 
the north-east direction from the main H Ii complex NGC 
2363 is seen in optical images. Since this paper treats 
the main star-forming region of a BCD, we concentrate 
on the main complex, which can be treated as a point 
source with the resolution of FIS. The north-east exten- 
sion is detected in the N60 (65 /.tm), WIDE-S (90 /im), 
and WIDE-L (140 /im) band images. This extended struc- 
ture does not affect our analysis because of its low surface 
brightness within the aperture used for the photometry. 

With IRAS data, Mazzarella, Bothun, & Boroson 
(1991) derived larger FIR fluxes (4.829 Jy at A 60 /im 
and 5.713 Jy at A = 100 //m) than our estimate. In the 
north-west extension, there is a brightness peak at dis- 
tance of ^ 3' from the main body. The flux around this 
peak is preliminarily estimated as 0.7 Jy at A = 90 /tm 
and 0.5 Jy at A = 65 /xm. This peak can be contaminated 
with the main body with the resolution of IRAS (~ 2'- 
5') (Neugebauer et al. 1984). Thus, we do not compare 
AKARI and IRAS results for this galaxy (section 3.2). A 
quantitative and detailed analysis of the north-east exten- 
sion will be reported in a separate paper. 

II Zw 70 — There is a bad pixel within the photometric 
radius in WIDE-S. Although the brightness of the pixel is 
extrapolated by using the value of the surrounding pixels, 
the uncertainty caused by this extrapolation may be as 
large as 30%.^ Thus, wc put a 30% error bar to the WIDE- 
S (90 /im) flux. 

II Zw 71 — This galaxy is occasionally detected in the 
same field of II Zw 70. It is located at the edge of the field 
of view, but the aperture used to determine the flux is well 
within the flcld. Although the area for the background is 
smaller than the other sources, the number of background 
pixels is large enough to subtract the background within 

^ We derived this error by using another analysis tool developed 
by Suzuki et al. (2008). 
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Fig. 1. Observed FIR SEDs at various bands (/^ is flux in units of Jy and A is wavelength in units of fira). AKARI data are shown 
by filled squares and IRAS data are also presented by open diamonds for eomparison. For II Zw 40, Spitzer data are available and 
are shown by open triangles. The bars indicate errors and the points with arrows give 3 a upper limits. The name of object is shown 
in each panel. 
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errors of 10% for N60 (65 ^JLu^) and WIDE-S (90 ^im), 
20% for WIDE-L (140 /Ltm), and 25% for N160 (160 ^m). 
This galaxy is not detected by IRAS. 

Mrk 36 — This is the faintest object in WIDE-S. This 
galaxy is detected by IRAS at A = 60 //m with a flux 
level of 0.022 ± 0.05 Jy and is not detected at A = 100 /ini 
with an upper limit of 0.68 Jy (Moshir et al. 1990). In 
the WIDE-L image, it seems that relatively bright pixels 
are clustered in the expected position of Mrk 36. If we 
use those clustered bright pixels and subtract the back- 
ground, the flux is 0.15 Jy. But since the uncertainty in 
the background is comparable to this value, it is safe to 
put an upper limit of 0.30 Jy. 

4. Derived Quantities 

4-.1. Dust temperature 

One of the most straightforward quantities derived by 
our observations is colors (i.e., flux ratios between two 
bands), which reflect dust temperatures. In Table 3, we 
show the colors. Since the WIDE-S (90 /^m) band is the 
most sensitive, we present the color relative to 90 ^m flux. 
We denote the flux ratio between wavelengths Ai [^m] and 
A2 [/^m] as [Al, A2]. The error of a flux ratio /1//2, where 
we denote the errors of /i and /2 as crr(/i) and crr(/2), 
respectively, is estimated by using a simple discussion of 
error propagation as 

We compare the derived colors with other available sam- 
ples of star-forming galaxies observed by AKARI FIS. To 
date, FIS observations of three spiral galaxies have been 
published to date (Kaneda et al. 2007; Suzuki et al. 2007) 
as shown in Table 3. We observe that the derived colors of 
the BCDs show higher temperatures (i.e., larger [65, 90] 
and smaller [140, 90]) than the spiral galaxies except for 
Mrk 7 and II Zw 71. This supports the idea that BCDs 
host intense star formation activities in concentrated re- 
gions. 

We have selected the present sample primarily from the 
IRAS catalog. Since intense star formation tends to raise 
FIR luminosity, our sample may be biased to objects with 
intense star formation. Indeed, II Zw 71, which is not 
detected by IRAS, has the lowest dust temperature, and 
this fact implies a possibility that our selection by IRAS 
detection is biased to high dust temperature. 

Now we estimate the dust temperature from the flux 
ratio. We fit Aiyl^B^{Td) (A is a constant, B^{T) is the 
Planck function evaluated at frequency v and tempera- 
ture T and P is the index of the emission coefficient called 
emissivity index) to the fluxes at two bands to derive the 
dust temperature Td. Note that A can be eliminated if we 
take a flux ratio. In Table 3, we show the estimated dust 
temperature, where Td(C; /3) represents the dust temper- 
ature derived from color C and emissivity index /3. 

From Table 3, we observe that the values of Td is higher 
for /3 = 1 than for /3 = 2, because under fixed fluxes at two 



wavelengths /3 = 2 requires lower Td (i.e., enhanced B^{Td) 
at the longer wavelength) to compensate for lower emis- 
sion coefficient at the longer wavelength. Again wc con- 
firm that the dust temperature of the BCDs is significantly 
higher than that of the spiral sample except for Mrk 7 and 
II Zw 71. We also see that there is a large variety of dust 
temperatures among our sample. Such a large variety of 
dust temperatures is also reported by Engelbracht et al. 
(2008) for a low-nietallicity sample. One of their sam- 
ple, II Zw 40 is common with Engelbracht et al. (2008), 
and the dust temperature derived by them by using the 
Spitzer 70 and 160 ^lm bands with /? 2 is 33.4 ± 1.0 K, 
consistent with our Td ( [140/90] ; = 2) = 36^^ K. 

If Td( [65/90]) is equal to rd( [140/90] ),6 the fluxes in 
the three bands can be fitted with a single dust tempera- 
ture. The difference between Td( [65/90]) and rd( [140/90]) 
is significant for Mrk 7 and II Zw 71 if /3 = 1, and for Mrk 7, 
Mrk 71, UM 439, UM 533, and II Zw 71 if /3 = 2. In those 
BCDs, Td( [65/90]) is higher than Td( [140/90]), which pos- 
sibly indicates that the contribution from stochastically 
heated grains to the 65 /Ltm flux is prominent compared 
with the short-wavelength tail of the emission at A > 90 fim 
(Draine & Li 2001). Since the emission at A > 90 iim 
comes from dust grains in temperature equilibrium be- 
tween the heating from stellar radiation field (called inter- 
stellar radiation field) and the cooling by FIR reemission, 
observations by the WIDE-S and WIDE-L bands are es- 
sential to estimate the equilibrium grain temperature. In 
such a case, the dust temperature derived from the IRAS 
60 fj,m and 100 fj,m bands is overestimate for the equilib- 
rium grain temperature. For BCDs in which Td( [65/90]) 
is consistent with Td( [140/90]), on the contrary, the IRAS 
60 fim and 100 fim bands have given reasonable estimate 
for the equilibrium dust temperature. In this case, the 
Spitzer 70 fim and 160 fim bands also provide the equi- 
librium grain temperature. 

It is worth noting that as demonstrated above the 
AKARI FIS bands have the following two advantages: 
One is that two bands, WIDE-S (90 ^im) and WIDE- 
L (140 /im), are available to estimate the equilibrium 
grain temperature, while the Spitzer 70 ^m and 160 
bands have risk of overestimating the equilibrium grain 
temperature because of the contribution from stochasti- 
cally heated grains to A = 70 /im. (But a wide band has 
a disadvantage in determining the flux density. In fact 
the color correction factor is relatively large for WIDE-S 
(90 ^m) and WIDE-L (140 //m) as stated in section 2.) 
The other is that the N60 band (65 /im) is also crucial 
to estimate the contribution from stochastically heated 
grains. Thus, the three bands of FIS adopted in this 
paper provide a unique opportunity to investigate both 
equilibrium-temperature dust and stochastically heated 
dust. 

^ When the parameter /3 is omitted, the statement is true for any 
0- 
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Name 


[65/90] 


[140/90] 


Td( [65/90] ;1) 


Td([140/90];1) 


rd([65/90];2) 


Td( [140/90]; 2) 


II Zw 40 


1.0±0.15 


0.56 ±0.12 


50^6 K 


53t'f K 


39I3 K 


3611 K 


Mrk 7 


0.99 ±0.29 


1.3 ±0.3 


48+1? K 


28+2 K 


38t^ K 


23tl K 


Mrk 71 


1.2±0.2 


0.71 ±0.15 


57+^8 K 


42+1" ^ 


43t^ K 


3lt^ K 


UM 439 


1.1±0.15 


0.78 ±0.17 


501^ K 


39ti" K 


39^3 K 


30t^ K 


UM 533 


0.89 ±0.12 


0.89 ±0.20 


44+^ K 


36t^ K 


35+^ K 


28l3 K 


II Zw 70 


1.1±0.35 


0.94 ±0.34 


53tll K 


34t^^ K 


4lt^° K 


27l^ K 


II Zw 71 


0.59 ±0.08 


2.4±0.5 


34+^ K 


22tl K 


29^2 K 




Mrk 36 


1.3±0.2 


0.58 ±0.58 


6ltl' K 


>30 K 




> 24 K 


NGC 2841 


0.29 ±0.08 


1.9 ±0.7 


24+^ K 


24+;^ K 


21+^ K 


20tt K 


NGC 2976 


0.43 ±0.12 


1.14±0.41 




K 


25ti^ K 


25+!^ K 


M 101 


0.36±0.10 


2.1 ±0.8 


27i^ K 


23t^ K 


23tl K 


19+t K 



Note. Lower limits are shown by > . ilie last three galaxies are spiral galax 
(2007) for NGC 2841 and NGC 2976 and from Suzuki et al. (2007) for M 101. 



4-2. Interstellar radiation field 

Since the dust temperature is determined by the equilib- 
rium between the heating from interstellar radiation field 
and the cooling by FIR reemission, it is possible to esti- 
mate the interstellar radiation field from the equilibrium 
dust temperature. We adopt Td( [140/90]) as the equilib- 
rium grain temperature. The 140 /im-100 ^m color of 
the Galactic disk is concentrated around 2.0 (Hibi et al. 
2006; Hirashita, Hibi, & Shibai 2007), which corresponds 
to the dust temperature of 20.2 K for (3 =\ and 17.3 K 
for (3 = 2. Assuming that the dust heating is dominated 
by UV radiation (Buat & Xu 1996), we estimate the UV 
radiation field from the dust temperature. The UV radi- 
ation field normalized to the Galactic value is denoted as 
Guv ■ Using the scaling relation of Guv oc '^t^^ (Evans 
1995), we obtain 



Guv(/3) 



Id 



20.2 K 

17.3 K 



if/3=l. 



if /? = 2. 



(2) 



In Table 4, we present Guv(/9) for P — 1 and 2. 
Naturally we confirm the above conclusion that the UV in- 
terstellar radiation field in BCDs is systematically higher 
than that in spiral galaxies, although Mrk 7 and II Zw 71 
are exceptions. Moreover, there is little difference in Guv 
between /3 = 1 and 2. Thus, it is more useful to adopt 
Guv instead of Td to avoid the dependence on (3. 

4-3. Total FIR luminosity 

Total FIR luminosity reflects how much stellar light is 
absorbed and reemitted by dust grains. It is empirically 
known that FIR luminosity is a good indicator of SFR 
in galaxies (Kennicutt 1998; Inoue, Hirashita, & Kamaya 
2000) . Here wc estimate the total FIR luminosity of each 
BCD. The FIS N60 (65 ^iui), WIDE-S (140 /im), and 
WIDE-L (140 /im) bands have the advantage of covering 
continuously the wavelength range from 50 /xm to 170 /im 
(Kawada et al. 2007). Thus, by using those three bands. 



we define the AKARI FIR flux, /akari, as an estimate of 
the flux in the above wavelength range: 



/akari = /.(65 ^iui)l\u{N60) 

+/^(90 ^lm)l:\v{WIDE-S) 
+/^(140 ^m)Ai^( WIDE-L) , 



(3) 



where /iy(A) is the flux per unit frequency at wavelength 
A, and Ai^(band) denotes the frequency width covered 
by the band. According to Kawada et al. (2007), we 
adopt Ai^{N60) = 1.58 THz, Av{WIDE-S) = 1.47 THz, 
and Av{WIDE-L) = 0.831 THz. With /akari, we define 
the AKARI FIR luminosity, Lakari, as an estimate of 
the luminosity in A ~ 50-170 fim: 



-^AKARi = 47r_D^ /akari , 



(4) 



where D is the distance to the object. In Table 4, we list 
the estimated AKARI FIR luminosities as well as the dis- 
tances. The distances are taken from Tolstoy et al. (1995) 
for Mrk 71 (NGC 2366), Osman, Ella, & Issa (1982) for 
NGC 2841, Karachentsev et al. (2002) for NGC 2976, and 
Jurcevic & Butcher (2006) for M 101. For more distant 
galaxies, we estimate the distance from the Galactocentric 
velocity (taken from NED)*" by assuming a Hubble con- 
stant of 75 km s^i Mpc^i. 

It is also possible to estimate the FIR luminosity in- 
cluding A > 170 fj,m by assuming a functional form of 
cx lyf'B^iTd), where we adopt Td = rd([140/90]). Then 
we obtain an estimate of the total FIR luminosity at 
A > 50 /im. Following the procedure in Nagata et al. 
(2002), we estimate the total FIR luminosity as below. 

First, we estimate the total FIR fiux. The total FIR 
fiux, /fir, can be well approximated by (Nagata et al. 
2002) 



/fir : 



^Tgo ( ) Bi,{Td)dv 



/ix(65 /im) - 5^(65 yum, Td) 



http : / / nedwww . ipac . caltech.edu / . 
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Table 4. Summary of the properties. 



Name D Guv -^akari -^fir 

(/3 = 1) (/3 = 2) (/3 = 2) (/3=1) (/3 = 2) 

[Mpc] [Lq] [Lq] [Lq] [Mq] [Mq] 

II Zw 40 9^2 1201^^" 8ll4^" 6.2 x 10^ 1.1 x 10^ 8.1 x 10*^ 3.9 x 10^ 1.1 x 10^ 

Mrk 7 42.4 5.4+2;o 5.7+2;o 1.9 x 10^ 2.3x10^ 2.2x10^ 1.4 x 10^ 3.7 x 10^ 

Mrk 71 3.4 39t^^ 34l^| 4.4x10^ 5.7x10^ 5.0x10^ 5.4 x 10^ 1.5 x 10^ 

UM 439 13.1 28t^4 251^^ 8.1x10^ 1.0 x 10^ 9.1x10^ 1.4 x 10* 3.7 x 10^ 

UM 533 10.9 17lf 16lf 8.1x10^ 9.9x10'' 9.1 x lO'' 2.2 x 10^ 5.8 x 10^ 

II Zw 70 17.0 141™ 14lf 2.9x10^ 3.4x10® 3.2x10® 8.1x10'' 2.1x10^ 

II Zw 71 18.5 1.4l:o4 1-5^0 5 2.5 x 10® 3.9x10® 3.5x10® 1.0 x 10^ 2.5 x 10^ 

Mrk 36 8.2 > 7.6 > 7.9 2.2 x lO'' — * — * > 8.2 x 10^ > 1.5 x 10^ 

NGC 2841 83 2.3tt;5 2^5^i^2 1.5 x 10^ 2.1 x 10^ 1.9 x 10^ 3.7 x 10^ 9.3 x 10^ 

NGC 2976 3.6 T.S+fg 8.1+^3 4.7 x 10® 5.8 x 10® 5.5 x 10® 3.1 x 10^ 8.1 x 10^ 

M 101 V'-^ l-Stgs 2.0tg;9 1.1 X 10™ .^-^^^0^°. 1-5x10^" 3.4 x lO'' 8 5 x lO'' 

Note. The last three galaxies are spiral galaxies for comparison (same as Table 3). Lower limits are shown by ">". 
* For Mrk 36, it is not possible to estimate LpiR because [140/90] is not determined. 



Xi2rgn 



-0 



AiyiN60), 



(5) 



where ft is the sohd angle of the object, tqq is the optical 
depth at A = 90 fiui, i/go is the frequency corresponding to 
A = 90 /xm (i.e., 3.33 THz), and B^(A, Tj) is the Planck 
function evaluated at wavelength A. The first term repre- 
sents the total FIR flux from grains in radiative equilib- 
rium, and the second term shows the contribution of the 
excess emission of stochastically heated grains in N60 (65 
/Ltm) . Thus, /fir provides the total FIR flux for A > 50 /xm. 
Noting that f^{90 fj,m) = i}TgoB^{90 fim, Td), we obtain 
the following expression similar to that of Nagata et al. 
(2002): 



/fir ^ Cif^{65 fim) 



/.(90 fim) 



exp 



159.^ 



- 1 



(6) 



exp(221.3/rd)-l_ 

where Ci = 1.58 x lO^^ Hz, Ca = 795.6 Hz K'^ and 
24.42 Hz K^^ for /3 = 1 and 2, respectively, and C3 = 
8.262 X 10^2 Hz and 1.144 x 10^^ Hz for /3 = 1 and 2, re- 
spectively. 

By using /fir, we obtain the total FIR luminosity ipiR, 

as 



-^FiR = 47r£> /fir , 



(7) 



where D is the distance to the object. In Table 4, we 
present the estimated FIR luminosities for the sample. 
The FIR luminosities estimated with (3 = 1 and 2 are in- 
dicated by Lfir(/3 = 1) and Lfir(/? = 2), respectively. For 
Mrk 36, it is not possible to determine ipiR since Td can- 
not be determined. Thus, we exclude Mrk 36 from the 
analysis here. The error of Lfir mainly comes from the 
uncertainty in the flux estimation (~ 20%). We note that 
the uncertainty caused by the difference in /? is compara- 
ble to or smaller than 20%. 



Table 5. Correction factor (Cakari) multiplied to the 
AKARI FIR luminosity to obtain the total FIR luminosity. 



(3 



1 
2 

fJote 



Spiral sample 
included not included 



1.36±0.17 1.34±0.19 
119 ±0 10 1.22 ±0.11 

'i'he number alter "±" shows the standard 



deviation (1 a). 



It is convenient to obtain a conversion formula from 
the AKARI FIR luminosity to the total FIR luminosity. 
We define the "correction" factor, Cakari, which is to be 
multiplied to the AKARI FIR luminosity to obtain the 
total FIR luminosity as 



LpmiP) = Cakari (/3)Takari , 



(8) 



where LpiR and thus Cakari depend on the assumed 
value of (3. Because the sample size is small, we cannot 
investigate the possible dependence of Cakari on various 
physical parameters such as dust temperature. Thus, we 
only estimate the average of Cakari for the sample. The 
results are listed in Table 5, where the averages of Cakari 
with and without the spiral sample are shown. We see 
that there is no significant difference in Cakari between 
those two averages. On the other hand, there seems to 
be difference in Cakari between /3 = 1 and 2, mainly be- 
cause the extrapolation to long wavelengths (A > 170 /im) 
in estimating the total FIR luminosity depends on (3. In 
summary, we can convert the AKARI FIR luminosity to 
the total FIR luminosity by multiplying 1.2-1.4. We note 
that a further 30-40% correction for the contribution from 
mid-infrared emission at A < 50 /im will be necessary to es- 
timate the total infrared luminosity of dust emission (Dale 
& Helou 2002). 
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4.4- Star formation rate 

It is empirically known that the FIR luminosity is a 
good indicator of the SFR. Hirashita. Buat, & Inouc 
(2003) provide a conversion formula from LpiR to SFR 
(we assume that Ljr in Hirashita. Buat, & Inouc 2003 is 
approximated by Lfir, but as mentioned in section 4.3, 
the upward 30-40% correction for A < 50 ^m may be re- 
quired for Lpm): 



SFR — CiR-LpiR , 



(9) 



where we adopt the conversion coefficient Cir = 2.0 x 
10~^° Mq yr~^ Lq^ as recommended by Hirashita, Buat, 
& Inoue (2003). The SFR estimated by using Lfir(/3= 1) 
is listed in Table 6 for each galaxy. If we adopt Lfir(/3 = 
2), we obtain smaller SFR by 4-26%. The range of the 
SFR is 0.01-0.4 Mq yr'^ for the present BCD sample. 
This is within the range (a few times 10~^ to several times 
10 Mq yr~^) given for a BCD sample by Hopkins, Schulte- 
Ladbeck, & Drozdovsky (2002). 

There are two common BCDs between this paper and 
Hopkins, Schulte-Ladbeck, & Drozdovsky (2002): Mrk 71 
and Mrk 36. For the latter galaxy, we do not obtain the 
SFR because it is not possible to determine Lfir (sec- 
tion 4.3). For Mrk 71, our estimate is 2-3 times lower 
than their estimate (0.065 Mq yr"^ from the IRAS 60 
fim luminosity and 0.041 Mq yr~^ from the 1.4 GHz ra- 
dio luminosity) . This confirms the comment of Hirashita, 
Buat, & Inoue (2003) that there is a risk of underestimat- 
ing the SFR for SFR < 1 Mq yr"! if we use equation (9). 
This may be because of smaller dust optical depth. 

There are three overlapping BCDs between this pa- 
per and Sage et al. (1992): II Zw 40, UM 439, and 
UM 533. Their SFR derived from the FIR luminosity 
(0.46, 0.085, and 0.049 Mq yi'^, respectively) is sys- 
tematically higher, which is attributed to different Cir 
(6.5 X 10~^ Mq yr"^ Lq'^). This high Cir is consistent 
with the SFR derived from Ha luminosity in dwarf galax- 
ies (Thronson & Telesco 1986), which again confirms the 
necessity of using higher Cir than proposed by Hirashita, 
Buat, & Inoue (2003) for SFR < 1 Mq yi'^ . Moreover, 
the SFR derived from Ha in Sage et al. (1992) (3.1, 0.20, 
and 0.19 Mq yr~^, respectively) is much higher than SFR 
derived from the FIR luminosity, which is again explained 
by small dust optical depth as mentioned in the previous 
paragraph. However, the SFR derived from Ha luminosity 
in Sage et al. (1992) could be overestimated because they 
assume that 2/3 of the ionizing photons are absorbed by 
dust, which might be less in low-metallicity galaxies like 
BCDs. 

4-5. Dust mass 

As stated in section 1, the WIDE-S (90 yum) and WIDE- 
L (140 fim) bands are suitable to trace the total mass of 
large grains which have the dominant contribution to the 
total dust content. The total dust mass Md is related to 
the flux as (e.g., Hildebrand 1983) 



M. = 



where is the mass absorption coefficient of dust grains. 
For simplicity, we assume spherical grains with a single 
radius a and with a uniform density s. As shown later, 
the assumption on a is not essential, since k^, is insensitive 
to a. The cross section of a grain is expressed as na^Q,^, 
where Qi, is a dimensionless quantity expressing the emis- 
sion efficiency. Using Qi/, the mass absorption coefficient 
is written as 

3Q, 



4as 



(11) 



Since there is a relation of Qi, oc a for a ^ A (Evans 1995), 
the estimate of is insensitive to the assumed value of 
a. Hildebrand (1983) estimates Q^/a at A = 125 jum to 
be 3/400 /dm~^. We note that the normalization of Qi,/a 
at A = 125 ^m is convenient since it is roughly the middle 
of the two wavelengths (A = 90 /.im and 140 fim) used for 
the estimate of dust mass. If we adopt s = 3 g cm^'^, we 
obtain 



18. 



A 



125 iim 



9 — 1 

cm g 



(12) 



(10) 



In Table 4, we list the dust mass, where we adopt 
Td([140/90]) for the dust temperature. The dust mass 
estimated with (3 = 1 and 2 are denoted as Md(/3 =1) 
and A/d(/3 = 2), respectively. The dust mass in Mrk 36 
is given by using the rd( [65/90]) instead of rd( [140/90]), 
for which only a lower limit is given. Since Td( [65/90]) 
tends to overestimate the equilibrium grain temperature 
(section 4.1), we regard this dust mass to be a lower limit. 
From Table 4, we observe that Md(/3 = 1) is systematically 
lower than Md(/3 = 2) by a factor of 2.5. This is because 
higher dust temperature for f3 = I leads to a lower dust 
mass with a fixed dust emissivity. 

One of our BCDs, II Zw 40, overlaps with the sample 
in Galliano et al. (2005) (section 3.3), who model the dust 
SED. First of all, the mass of the big grain component, 
which is contributing to the FIR regime, is just the same 
as our estimate (1.1 x 10^ Mq) if we adopt (3 = 2. This 
means that the AKARI WIDE-S (90 /im) and WIDE- 
L (140 yum) bands are indeed useful to trace the dust 
mass. However, they also show an additional very cold 
dust component contributing to the emission at submil- 
limeter wavelengths, although the amount of this compo- 
nent is dependent on the FIR emissivity index /3. This 
very cold component cannot be investigated by AKARI, 
and we concentrate on the dust component contributing 
to the FIR emission in this paper. 

5. Discussion 

5.1. Star formation properties 

In section 4.2, we have shown that the BCD sample 
tends to show higher UV radiation fields than spiral galax- 
ies. This is consistent with a picture that intense star for- 
mation is occurring in a concentrated region in BCDs. In 
order to examine the properties of such an intense star 
formation activity (e.g., a burst or a continuous mode of 
star formation) it is useful to estimate the gas consump- 
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Table 6. Some related quantities used for discussion. 



Name 


or xl 


iWHi 


''SF 


19^ Inn- (Ct /TT^ 


T) 




[Mq yr J 










TT 7ttt /in 


U.zz 


z.U X iU 


u.yi 


o.io 


z.U X lU 


Mrk 7 


0.50 


3.6 X 10^ 


7.2 


8.54 


3.9 X 10-'' 


Mrk 71 


0.012 


1.2 X 10^ 


110 


7.83 


4.5 X 10"^ 


UM 439 


0.021 


1.7 X 10^ 


8.1 


7.98 


8.2 X 10-'"^ 


UM 533 


0.020 


5.8 X 10^ 


2.9 


8.10 


3.8 X 10-4 


II Zw 70 


0.073 


3.6 X 10^ 


4.9 


8.11 


2.3 X 10-" 


II Zw 71 


0.080 


7.3 X 10^ 


9.1 


8.24 


1.4 X 10-3 


Mrk 36 


_t 


1.8 X 10^ 


<2.5 


7.81 


>4.6x 10-5 



Note. Lower and upper limits are shown by ">" and "< ", respectively. Vox Mrk yti, bl<'R 
cannot be estimated. 

* The SFR is derived by using -Lfir(/3 = !)■ If ^firX/S = 2) is used instead, the SFR becomes 
4—26% higher according to the difference in Lfir. 

t The SFR of Mrk 36 cannot be estimated because Lfir is not available. 



tion timescale (tsf) as follows. 

Here, we adopt the SFR estimated in section 4.4 by 
adopting Lfir(/3 = 1). (If we adopt Lfir(/3 = 2), the gas 
consumption timescale becomes 4-26% shorter.) Then 
the gas consumption timescale is defined as the gas mass 
divided by the SFR. For the gas mass, we adopt the 
mass of Hi gas (A/jji) in the same way as in Hirashita, 
Tajiri, & Kamaya (2002) and Hopkins, Schulte-Ladbeck, 
& Drozdovsky (2002), since molecular gas is rarely de- 
tected for BCDs (Barone et al. 2000). We estimate Mri 
of II Zw 71 by using the Hi flux measured by Huchtmeier 
& Richter (1989) and by using equation (1) in Lisenfeld 
& Ferrara (1998). The estimated gas consumption time 
TsF is listed in Table 6. We observe that tsf is dis- 
tributed with a range of two orders of magnitude, which is 
consistent with Hopkins, Schulte-Ladbeck, & Drozdovsky 
(2002). Moreover, the gas consumption time does not 
correlate with Guv- For example, Mrk 71 has large Guv, 
which indicates a concentrated intense star formation, but 
has large tsf exceeding the cosmic timescale. Such galax- 
ies with long Tsf but high Guv give us a picture that the 
star formation is strongly activated only in a limited re- 
gion. On the contrary, II Zw 40 has large Guv and small 
Tsf, indicating the whole galaxy may be in an active phase 
of star formation. 

It is possible that the star formation in BCDs is in- 
termittent. In this case, the diversity in tsf is naturally 
explained: if a BCD is in activated/inactivated phase of 
star formation, tsf becomes small/large. Indeed, the ef- 
fect of feedback is considered to be enhanced in small sys- 
tems such as BCDs, leading to an intermittent nature of 
star formation (Saito, Kamaya, & Tomita 2000; Carraro 
et al. 2001; Kobayashi & Kamaya 2004; Kamaya 2005). It 
is also observationally suggested that the star formation 
in dwarf galaxies may be episodic (Fanelli, O'Connell, & 
Thuan 1988; Greggio et al. 1998; Grossi et al. 2007). 

5.2. Dust-to-gas ratio and chemical enrichment 

The dust production is strongly related to the metal 
production. Indeed, the relation between dust abundance 
(usually dust-to-gas ratio is adopted) and metal abun- 



dance (metallicity) can be investigated by using chem- 
ical evolution models (Dwek 1998; Hirashita, Tajiri, & 
Kamaya 2002). Here we consider whether dust-to-gas ra- 
tio and metallicity are related or not. 

As the metal abundance, we adopt oxygen abundance, 
which is the most easily observed. The oxygen abun- 
dance of each sample BCD is compiled in Hirashita, 
Tajiri, & Kamaya (2002) and Hopkins, Schulte-Ladbeck, 
& Drozdovsky (2002), but for II Zw 71 we take the oxygen 
abundance from Kewley, Jansen, & Geller (2005). The 
dust-to-gas ratio T> is estimated by M^iP = 1) divided by 
Mb_i. If we use M^i/S = 2) instead, we obtain about 2.5 
times higher dust-to-gas ratio. Note that the dust-to-gas 
ratio would shift further upward if there is really a very 
cold dust component contributing to longer wavelengths 
than those traced by AKARI (Galliano et al. 2005; sec- 
tion 4.5). In Table 6, we list the metallicity (oxygen abun- 
dance: 12 -I- log(0/H)) and the dust-to-gas ratio V. Note 
that the solar abundance is 12-|-log(0/H) = 8.93 (Anders 
& Grevesse 1989). 

In Figure 2, we show the relation between those two 
quantities. We observe that there is a positive correlation 
between dust-to-gas ratio and metallicity. This correlation 
is also shown by Engelbracht et al. (2008) with a larger 
metallicity range. The correlation is consistent with the 
picture that dust enrichment proceeds as the system is 
enriched by metals. 

It is also interesting to compare the results with the- 
oretical model predictions. Here we adopt a model by 
Hirashita, Tajiri, & Kamaya (2002) (originally developed 
by Lisenfeld & Ferrara 1998), who consider a one-zone 
chemical evolution model with the instantaneous recycling 
approximation (Tinsley 1980). The model equations con- 
sist of the evolutions of gas mass, metal mass, and dust 
mass. By treating these equations, the relation between 
dust-to-gas ratio and metallicity is obtained. In partic- 
ular, there are two important parameters in the model: 
(i) the fraction of metals condensed into dust grains, /in, 
and (ii) the "dust destruction efficiency", Ps'n- The for- 
mer parameter is fixed to 0.1, and does not contribute 
to the scatter of the dust-to-gas ratio if the condensa- 
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Fig. 2. Relation between dust-to-gas ratio "D and metallic- 
ity. The data of the present BCD sample are shown by solid 
squares. The arrow indicates a lower limit. The solid, dotted, 
and dashed lines show the model predictions by Hirashita, 
Tajiri, & Kamaya (2002) with various dust destruction effi- 
ciencies (in their notation /3sN = li 5, and 25, respectively). 

tion efficiency in the stellar mass loss is independent of 
age. The "dust destruction efficiency" /3sn is defined as 
PsN = Mg/{TsNij), where Mg is the total gas mass, tsn 
is the timescale of dust destruction by supernova shocks, 
and tjj is the SFR. This parameter can be expressed as 
follows (Hirashita, Tajiri, & Kamaya 2002, equation 6): 



7 1 
/3sN = eAfs — - 



where e is the fraction of destroyed dust in a SN blast, 
Ms is the gas mass accelerated to a velocity large enough 
for dust destruction (typically 100 km s~^), 7 is the su- 
pernova rate, and Xgp is the gas mass fraction in the 
star-forming region, where the dust formation and de- 
struction are really occurring (the rest of the gas is con- 
sidered to be contained in the H i envelope, which they 
assume not to be concerned with the dust formation and 
destruction). Typical values are e ^ 0.1 (McKce 1989), 
Ms ~ 68OOM0 (if we adopt 10^"'^ erg for the kinetic energy 
input of a supernova, and 100 km s^^ for the threshold 
velocity for the dust destruction), j/ip 1/136Mq (for 
the Galaxy; Lisenfeld & Ferrara 1998), we obtain /3sn = 5 
for XsF = 1- Hirashita, Tajiri, & Kamaya (2002) argue 
that this parameter changes according to the relative con- 
tribution between Type la supernovae and Type II su- 
pernovac. If an intermittent star formation activity is as- 
sumed, j/i/j can change by a factor of 20 (Bradamante, 
Matteucci, & D'Ercole 1998; Hirashita, Tajiri, & Kamaya 
2002). Motivated by this factor, they changed /9sn from 
1 to 25. 

In Figure 2, we show the theoretical predictions for 
/? = 1, 5, and 25 with /in = 0.1. We observe that those 
models explains the data points except for the two ex- 
tremes (II Zw 71 and Mrk 71 for the upper and lower 



extremes in Figure 2, respectively). The variation in /3sn 
increases the scatter of the dust-to-gas ratio at a certain 
metallicity and makes the correlation less clear. Further 
change of /?sn may be possible because of the change of 
XsF- As mentioned in section 5.1, only a small fraction 
of gas may be concerned with star formation in Mrk 71. 
Thus, XsF <SC 1 is expected for Mrk 71, where /9sn could 
be larger than the above values. Such a large value of /?sn 
may explain the extremely poor dust-to-gas ratio of Mrk 
71. 

As commented above, the dust-to-gas ratios are larger 
if we adopt (3 = 2, and/or if the cold dust component 
contributing to the submillimeter wavelengths, which we 
fail to trace in FIR, really exists as proposed by Galliano 
et al. (2005). In these cases, we should adopt a larger /in, 
but a large variation of /3sn is still required to explain the 
large scatter of dust-to-gas ratios. 

6. Summary 

We have reported far-infrared (FIR) properties of eight 
blue compact dwarf galaxies (BCDs) observed by AKARI. 
The fluxes at wavelengths of 65 /im, 90 jum, 140 jum, 
and 160 /im are measured and arc used to estimate ba- 
sic quantities about dust grains such as dust temperature, 
dust mass, and total FIR luminosity. We have found that 
the typical dust temperature of our BCD sample, except 
for Mrk 7 and II Zw 71, is significantly higher than that 
of normal spiral galaxies. This indicates high interstel- 
lar radiation field, which proves to range up to 100 times 
the Galactic value. This confirms the concentrated star- 
forming activity in BCDs. The star formation rate (SFR) 
is also evaluated from the FIR luminosity as 0.01-0.5 Mq 
yr~^. Combining this quantity with gas mass taken from 
the literature, we estimate the gas consumption timescales 
(gas mass divided by SFR), which span a wide range from 
1 Gyr to 100 Gyr. A natural interpretation of this large 
variety can be provided by intermittent star formation ac- 
tivity. We finally show the relation between dust-to-gas 
ratio and metallicity (we utilize our estimate of dust mass, 
and take other necessary quantities from the literature). 
There is a positive correlation between dust-to-gas ratio 
and metallicity as expected from chemical evolution mod- 
els. 
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